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ABSTRACT 

We investigate the large-scale structure of the local galaxy distribution using the 
qq ■ recently completed 2 Micron All Sky Survey (2MASS) via three techniques. First, 

we determine the AVband number counts over the «4000 deg 2 APM survey area 
where evidence for a large-scale 'local hole' has previously been detected and compare 
them to a homogeneous prediction. Considering a ACDM form for the 2-point angular 
correlation function, the observed deficiency represents a 5cr fluctuation in the galaxy 
distribution. We check the model normalisation using faint /\-band data compiled from 
the literature; the normalisation used in this paper is in excellent agreement, and the 
observed counts over the APM survey area would require the model to be lowered by 
3.8<t. However, the issue is complicated by the b ^20° and b ^-20° 2MASS counts 
which lie below the best-fit model normalisation. Second, since the -ft^-band counts 
over the APM survey area continue to suggest the possible presence of excess clustering 
over the ACDM prediction, we next probe the power at large scales by comparing the 
2MASS and ACDM mock galaxy angular power spectra. We find a 3cr excess in the 
q ■ 2MASS catalogue over the ACDM prediction at large scales (I ^ 30). However, this 

?-h ' excess is not enough to account for the low counts over the APM survey area. Finally, 

, in order to probe more directly whether the ACDM mocks can reproduce observed 

features in the galaxy distribution at large scales, we apply a counts in cells analysis to 
the 2MASS data and mock catalogues; on the assumption that the 2MASS catalogue 
at \b\ ^ 20° is representative, we find excellent agreement between the biased ACDM 
mocks and the 2MASS catalogue to ~ 30°. The crux of the interpretation of these 
results appears to be whether the 2MASS volume is yet big enough to constitute a fair 
sample of the Universe. The number count models based on fainter counts suggest that 
it may not be, although normalisation uncertainties remain. It is also the case that the 
2MASS depth remains comparable to the possible size of large-scale inhomogeneities. 
Analyses which assume the 2MASS average density is a fair sample, such as counts in 
cells and to a lesser extent, power spectral analysis, may return results which see less 
contradiction with ACDM than the number count analysis in the APM survey area. 
Further progress on assessing the significance of the local hole and thus the consistency 
of the local galaxy distribution with ACDM will require deeper all-sky number counts 
and redshift surveys in both the optical and the infrared. 

Key words: surveys - galaxies: photometry - cosmology: observations - large-scale 
structure of the Universe - infrared: galaxies 
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1 INTRODUCTION 

The counting of galaxies as a function of apparent mag- 
nitude is one of the most powerful tools in observational 
cosmology. Not only can this simple statistic form strong 
constraints on the level of evolution at the faint end, but 
also on the large-scale structure and the scales to which the 
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cosmological principle can be said to hold from bright mag- 
nitude counts. 

The APM galaxy survey dMaddox et alll99tj|) produced 
number counts which are unexpectedly low for the solid an- 
gle surveyed («4000 deg 2 ). The question then arises: To 
what extent are these steep counts due to real clustering 
in the galaxy distribution, and do errors in the photomotery 
or strong evolution at low redshifts contribute significantly 
to the observed deficiency? If the counts were exclusively due 
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Driver! 120031) and the Sloan Digital Sky Survey (SDSS; 



Yasuda et all 1200 J) Early Data Release. However, a com 



Figure 1. Number redshift histograms for the 2MASS-2dFGRS 
matched sample described in section 2.2 for the Northern (left- 
hand panels) and Southern 2dFGRS declination strips (right- 
hand panels). The upper panels show the n(z) for K s <13.5 
galaxies. The lower panels show the correspond ing radial den- 
sity fu nctions, i.e. the n(z) divided through by the|Metcalfe^^^ 



2001) homogeneous prediction. In each case the Metcalfe ct al. 
2001) model is indicated by a solid line. We also show the la 



lower limit best fit normalisation of the Metcalfe ct al. (20Q1J) K- 
band model detailed in section 3.1 (dot-dashed lines). 



to an under-density in the Southern Galactic Cap (SGC), it 
would be unexpectedly large and deep for our present under- 
standing of large-scale structure. However, to invoke galaxy 
evolution alone is also problematic as it requires an extended 
tail in the n(z) which is not apparent. 

The presence of a large under-density in the SGC has 
been confirmed by large galaxy re dshift surveys. The 2dF 
Galaxy Redshift Survey (2dFGRS: IColless et al.ll200ll) has 
taken spectra for galaxies brightward of bj «19.45 over a 
solid angle of >600 deg 2 in the SGC. The n(z) indicates 
remarkable structure in the local galaxy distribution with 
a large deficit to z ~0.1. Large deficiencies in the n(z) 
are also indicated by other redshift surveys in the SGC 
jRatcliffe et al.lll998t IShectman et alJll996t IVettolani et alJ 
119971) . 

Possible errors in the APM pho tometry have al s o been 
detected. Using CCD photometry, iMetcalfe et al] dl995l) 
found a small residual scale error in the APM survey zero- 
points for B> 17. The corrected counts fainter than B ~17.5 
were now in good agreement with the expected homogeneous 
trend. However, a deficit in the counts of «50 per ce nt at 
B ~16 remained. More recently. iBusswell et al] l|2004l) have 
provided B-band CCD photometry over f»337 square de- 
grees in the SGC to a limiting magnitude of B=20.5. The 
photometric zero-point is in excellent agreement with other 
CCD data such as the Millennium Galaxy Catalogue (MGC; 



parison with the APM Bright Galaxy Catalogue suggests 
a large zero-point offset of 0.31 magnitudes for B < 17.35. 
Appl ying this correc tion to the bright APM photometry, 
IBusswell et al] (|2004) found that the APM survey counts 
now indicated a ~25 per cent deficit at B «16, consistent 
with a 25 per cent deficiency in the galaxy distribution to 
2=0.1 over the entire APM survey area. The implied sig- 
nificance using several possible forms for the real-space cor- 
relation function indicated that such a large local hole in 
the galaxy distribution required excess power at large scales 
over the ACDM prediction. 

In conjunct io n with the work of IBusswell et all ]2004) . 
iFrith et al] ]2003l) investigated the number counts in the 
AVband from the 2 Micr on All Sky Survey (2MASS) sec- 
ond incremental release Jjarrett et al] l2000l) . They found 
good agreement both with the optical number counts, and 
the expected trend defined by the 2dFGRS n(z), suggesting 
that the number counts in the 2dFGRS strips are caused 
by real features in the galaxy distribution, and are consis- 
tent with an absence of strong, low-redshift galaxy evolution. 
However, due to incompleteness in the 2MASS second incre- 
mental release, a comparison with the optical APM survey 
counts could not be made, although the steeper than Eu- 
clidean counts over the partially surveyed galactic caps in- 
dicated the possible presence of large under-densities in both 
the Southern and Northern Galactic Caps to ~300 ft. -1 Mpc. 
This was supported by mapped 2MASS counts, showing the 
variations of counts over the sky in 5°x5° bins, and indi- 
cating huge regions of underdensity > 100° across in both 
hemispheres. 

Here, we probe the large-scale structure of the local 
galxy distribution with the recently completed 2MASS sur- 
vey using three techniques. First, we compare the A's-band 
2MASS counts over the APM survey area with model pre- 
dictions. We calculate the significance with respect to these 
models considering a ACDM form for the 2-point angular 
correlation function . Secondly, since the form of the local 
hole as suggested bv lBusswell et aD (|2004l) requires the pres- 
ence of excess power at large scales over the ACDM pre- 
diction, we compute the 2MASS angular power spectrum 
and compare with a ACDM prediction constructed from the 
ACDM Hubble Volume simulation. Finally, we investigate 
the local galaxy distribution more directly by using a counts 
in cells analysis on the 2MASS and ACDM mock catalogues. 

In section 2, we present details of the datasets used. In 
section 3, the number counts are presented. The clustering 
at large scales is investigated through a determination of the 
2MASS and ACDM mock angular power spectra in section 
4. In section 5 we present the counts in cells analysis. The 
discussion and conclusions follow in section 6. 



2 DATA 

2.1 The 2MASS Extended Source Catalogue 

The 2 Micron All Sky Survey (2MASS) extended source 
catalogue has now released K s , H and J-band photometry 
for over 1.6xl0 6 extended sourc es over the entire s ky with 
high completeness to 7^=13.5 ]Jarrett et aljEood) . Since 
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the second incremental release, much of the photometry has 
been revised (Jarrett - priv. comm . ) and the default aper- 
ture m agnitudes u s ed inlCole et alJ (|200lt) have been aban- 
doned. ICole et aD i200ll) compared the second incremen- 
tal rel ease -RT s -band p hotometry with the TVband photome- 
try of lLovedavl (l200(tl . The Loveday photometry had better 
signal-to- noise and resolution than the 2MASS scans and 
so enabled more accurate 2MASS magnitudes to be deter- 
m ined. Here, we us e a similar magnitude estimator to that 
of lCole et alJ <l200ll) . We take the J- band extrapolated mag- 
nitudes, colour- corrected to the K 3 -band using the J and 
-ff s -band fiducial elliptical magnit udes; this magn itude esti- 
mator provide s the best fit to thelLovedavl feOOOT) zeropoint 
as described in IFrith et all d2004bl) An extinction correction 
is applied using the ISchlegeTetHd] l)l998l) dust maps. 

2.2 The 2dF Galaxy Redshift Survey 

The 2dF Galaxy Redshift Survey (2dFGRS) is selected in 
the photographic fej-band using the A PM survey and sub - 
sequent alterations and extensions to it dColless et alJl200lh 
for two declination strips in the northern and southern galac- 
tic caps, as well as 99 randomly selected 2° fields scattered 
over the APM survey area. The final release data has pro- 
vided spectra for ~220,000 galaxies (for a quality flag of 
Q^3) over ~1500 square degrees to a mean magnitude limit 
of 6j = 19.45. 

In this paper, we investigate the galaxy distribution in 
the .Ks-band. We have therefore formed a 2MASS-2dFGRS 
catalogue, matched over the Northern and Southern 2dF- 
GRS declination strips, using the K s -band magnitude esti- 
mator described above. Fig. shows the iG-band selected 
n(z) and radial density functions, i.e. the observed n(z) over 
the predicted, for the 2dFGRS declination strips with an ap- 
plied limiting magnitude of K s =13.5. 

2.3 The ACDM Hubble Volume Simulation 

The Hubble Volume catalogues represent the largest N-body 
simulations of the Universe to date. The ACDM simula- 
tion follows the evolution of 10 9 dark matter particles from 
z «50 over a volume of 3000 3 /i _3 Mpc 3 . The associated cos- 
mological parameter s are fi m =0.3, f2{,=0.04, h—0.7, as=0.9 
J Jenkins et al.lll998h . 

In this work, we construct mock 2MASS catalogues 
from the z = ACDM Hubble Volume simulation. We di- 
vide the total volume into 27 virtually independent spheri- 
cal volumes of r — 500 h^ 1 Mpc. These are subjected to the 
2MASS selection function: 

H{Z) = 2{zJlA12f eXP (- {^T) (1) 

l|Baugh fc Efstathiodll993l: iMaller et aLll2003ah where z is 
determined from the 2MASS-2dFGRS matched sample de- 
scribed in section 2.2. Equation is normalised to match 
the total number of observed 2MASS galaxies for \b\ ^20°. 
Due to the volume of the 27 mock 2MASS catalogues, the 
selection function is artificially truncated for the K s <13.5 
mocks at z ~0.156. However, this has a negligible effect on 
the work in this paper; at this redshift, ~95 per cent of the 
galaxies are sampled for K s <13.5. 



For the counts in cells analysis, it is necessary to more 
accurately mimic the galaxy sample. We therefore use a bias 
prescription: 

T>(,A - / ex P( a ^ + Pv 3/2 ) for v ^ , . 

1 ' ~ \ exp(c^) for v < 0, { > 

\Co\e et al.lll998l) where P{v) is a bias probability based on 
the density field and v is the number of standard deviations 
of the density away from the mean. For the counts in cells 
analysis we use magnitude-limited data of K s <12.5. We use 
parameters of a— 0.6 and /3= -0.15 to match the K s <12.5 
angular correlation function at small scales (6< 10°). 



3 NUMBER COUNTS 

3.1 Model Normalisation 

The issue of the local hole rests critically on the model 
number count normalisation. In the K-band, the num- 
ber count predictions for passively-evolving galaxies are 
in remarkable agr e ement with the observations to K «23 
iMcCracken et all |2000). The predicted number below 
K «18 is also fairly insensitive to the evolutionary model 
or the assumed cosmology. Therefore, the K-band num- 
ber counts are a particularly useful probe of the local Uni- 
verse since the model predictions can be constrained at 
fainter magnitudes with few concerns over uncertainties in 
the amount of evolution or the cosmology. 

In this work, we use a non-evolving TVband 
model computed from the l uminosity function parame- 
ters of iMetcalfe et al.| ll200ll) and the if-corrections of 
iBruzual fc Charlotl Jl993l) . Comparing the number count 
predictions in the fitt ing range 14< K <1 8 to the 
faint counts u sed in iMcCracken et al] j20 0Cf) as well 
as other data tVaisjmgjiiyLalJ 120001: ISzokolv et alJ Il998l; 
Huang et alJll997l. l200ll; iKummel fc Wagnei|l2000l iMartinil 
200(j iMcLeod et alJll995l) . we find that the lMetcalfe et alJ 
model provides a good fit to the ob servations with 
y=0.9 6±0.06 (where Y=l represents the IMetcalfe etafl 
i200ll) normalisation). 

In order to test for the presence of strong evolution at 
low redshifts and zero-point offsets in the 2MASS data, we 
als o construct model pred ictions for the number counts from 
the lMetcalfe et al.l J2001) hom ogeneous predictio n described 
above and the observed n{z) iFrith et alJl2003l) . These are 
determined by varying the luminosity function normalisa- 
tion as a function of the redshift; the luminosity function 
parameter <f>* is multiplied by the relative density (Fig. 
lower panels). These variable <fi* models therefore provide a 
simple prediction for the number counts for a given galaxy 
distribution, assuming that there is no significant effect from 
strong galaxy evolution at low redshifts. 

3.2 Results 

In order to verify the consis tency of bright n umber counts 
with the corresponding n{z) iFrith et al.ll2003T) . we compare 
the iCj-band number counts extracted from the 2dFGRS 
strips with the corresponding variable (f>* models, using the 
newly-completed 2dFGRS and 2MASS catalogues (Fig. 0. 
The .Ks-band number counts extracted from the 2dFGRS 
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Figure 2. its-band 2MASS number counts extracted for the 2dF- 
GRS Northern (left-hand panels) and Southern (right-hand pan- 
els) declination strips. The upper panels show the counts on an 
ordinary logarithmic plot, while the lower panels show the resid- 
ual plo t of the counts , i.e. the number count divided through 
by the iMetcalfe et alj 1200 J) homogeneous prediction. The ho- 
mogeneous model is indicated by a solid line, with the 2dFGRS 
Northern and Southern variable (f>* models shown by the dotted 
and dashed lines respectiv ely. We also show the la lower limit 
best fit normalisation of the lMetcalfe et all l200lf) it-band model 
detailed in section 3.1 (dot-dashed lines). The errorbars are Pois- 
sonian. 



fields are in reasonable agreement with the corresponding 
variable <j>* models. This suggests that the form of the bright 
number counts is exclusively due to features in the local 
galaxy distribution and is consistent with an absence of 
strong galaxy evolution at low redsbifts. The agreement be- 
tween the counts and the variable cj>* models is independent 
of the model normalisation, since any change in the number 
count models also alters the n(z) model normalisation and 
therefore the implied deficiency to the same degree. 

Having confirmed the consistency between the number 
counts and the underlying large-scale structure in the 2dF- 
GRS strips, we are now in a position to examine the number 
counts over the APM survey area. Fig. [3]sliows the it s -band 
2MASS counts extracted for the f«4000 deg 2 field. We also 
show the 2dFGRS Southern variable </>* model for reference. 

There is a surprisingly good agreement between the 
its-band counts and the 2dFGRS Southern variable <j>* 
model. This suggests that the local galaxy distribution in 
the APM survey area may be similar to that of the much 
smaller 2dFGRS Southern strip. This is supported by sim- 
ilar def iciencies observed in the Durham-UKST redshift 
survey (|P^tdiffe^t^ljJT998) , the Las Campanas Redshift 
Survey Jsiiectman^t^nTl99a) and the ESO Slice Project 
iVettolan^^dT^97l) ~whicIi are all situated within the 



Figure 3. it s -band 2MASS number counts extracted for the m 
4000 deg 2 APM survey area and the \b\ 20° galactic caps, 
presented as in Fig. [5] 



APM survey area, and also the optic al number counts 
and corresponding variab le <j>* models (iFrith et al . 2004a, 
l2003tlBusswell et all20ol . The degree of under-density may 
therefore be inferred from the observed 2MASS-2dFGRS 
redshift distribution and th e model normalisation. Tak- 
ing the lMetcalfe et ail teOOll) it-band model normalisation 
and the 2MASS-2dFGRS matched n(z) (Fig. [TJ implies a 
24 per cent under-density to z=0.1. If we take the la lower 
limit for the best fit it-band model normalisation detailed in 
section 3.1, the implied under-density over the APM survey 
area to z=0.1 is 15 per cent. 

In order to see whether this structure persists to even 
larger scales than the APM survey area, we have determined 
the 2MASS it s -band counts over the entire |6| ^20° galac- 
tic caps (Fig. |3J. The counts are very low with respect to 
the model normalisations used previously, but are in good 
agreement with ea ch other. This may ind icate that the nor- 
malisation of the IMetcalfe et alJ i200ll) model should be 
even lower, 2.5cr below the best fit to the it-band counts at 
14< it <18. Alternatively, it might also indicate the pres- 
ence of a zero-point offset between the 2MASS photometry 
and the it-band model of ~0.15 magnitudes; any alteration 
to the 2MASS zero-point would compromise the agreement 
between the it s -band 2dFGRS strip counts and the cor- 
respon ding variable (f>* models, and also with the iLovedavl 
(2000) photometry zero-point with which the 2MASS mag- 
nitudes are compared. If the low counts over the galac- 
tic caps were due to real features in the galaxy distribu- 
tion, this would imply that the local Universe is globally 
under-dense, and that even surveying over «25,000 deg 2 to 
r> 150 h~ x Mpc does not constitute a fair sample of the Uni- 
verse. 
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Model Normalisation 


Y 


Deficiency 


Significance 


Metcalfe et al. 


1.00 


26.9 per cent 


5.0a 


Best fit 


0.96 


23.9 per cent 


4.5(7 


1<t lower limit 


0.90 


18.8 per cent 


3.5(7 


" 2MASS 


0.80 


9.7 per cent 


1.8(7 



Table 1. Significances for the observed deficiency in the X 3 -band 
counts extracted for the APM survey area for K B <12.5 consid- 
ering various model normalisations (the relative normalisations 
are indicated by Y). In each case a ACDM correlation function is 
considered of A=0.28 l-7=-0.71 for 9 <5.0°, A=5.3, l-7=-2.5 for 
8 Jj5.0° where u j = AO 1 f . The best fit parameters at small scales 
are taken from iMaller et al ] i2003al) . while the large scale fit is 
determined from the K s <12.5 correlation function presented in 
Fig.Elfor 9 ^10°. This is formed from the 27 mock 2MASS cata- 
logues constructe d from the AC DM Hubble Volume mocks with 
a bias of b ft =1 .1 jMaller - e^ rt.ll200.Sa|) and is in good agreement 
with the lMaller et alj \2003e[ power law best fit at 6 ml° . 

3.3 Determining the Significance 

Using assumed forms to the galaxy correlation function at 
large scales, it is possible to determine the as sociated signif- 
icance of features in the galaxy distribution. iBusswell et alJ 
(2004) considered various forms to the real-space correla- 
tion function to calculate the significance of an assumed 
3-dimensional form to the local hole. Here, we use the 2- 
dimensional analogue to determine the significance implied 
by the observed deficiency in the n umber counts w i th re - 
spect to various normalisations of the lMetcalfe et ail (1200 if) 
model using the angular correlation function: 

(f ) 2 = i + ^/^ 2 -(M (3) 

JPeeblesl ll980l where uj (612) is the value of the two-point 
angular correlation function between two area elements dQ,± 
and dQ2- f2 is the total solid angle of the survey, n is the 
mean galaxy density such that N = nil is the total num- 
ber of galaxies in the survey area. A power law form of the 
correlation function is assumed such that u) = AO 1 ^ 1 . 

In Table we apply this technique to the observed 
deficiency over the APM survey area in t he 2MASS K s - 
band counts for K s <12.5 with respect to the lMetcalfe et ail 
(2001) model, the best fit and the lcr lower limit best fit 
normalisations of this model described in section 3.1, and 
the observed mean 2MASS number count for \b\ ^ 20° (see 
Fig. [3 • This magnitude limit is chosen such that the peak 
in the selection function lies within the redshift range of in- 
terest; at K a <12.5 the n(z) peaks at z «0.05, and so this 
magnitude limit preferentially samples the mean depth of 
the possible local hole in the APM survey area. We consider 
a ACDM form for the angular correlation function deter- 
mined from fits to the mock 2MASS catalogues (see section 
4.3). 

Both the lMetcalfe et alJ feOOlfl model and the best fit 
normalisation appear to imply that the observed 2MASS 
number counts over the APM survey area represent either 
an extremely rare fluctuation in the galaxy density or a chal- 
lenge to the large-scale clustering predicted by the ACDM 
Hubble Volume simulation. Even lowering the normalisa- 
tion to the lcr b est fit lower limit yie lds a significance of 
> 3cr. Only if the lMetcalfe et all feOOlft model is lowered by 



2.5a to the mean 2MASS number density at K s <12.5 for 
\b\ ^ 20° (n 2mass) , do the counts become more consistent 
with the ACDM prediction. To remove any discrepancy in 
the K a <12.5 counts over the APM survey area requires the 
model normalisation to be lowered by 3.8cr. Combining the 
error on the normalisation and the significance estimate for 
the best fit normalisation in quadrature yields a significance 
of 3.0cr. Clearly, unless the faint i\s-band counts are sys- 
tematically too high, then the 2MASS counts over the APM 
survey area suggest the presence of excess power at large 
scales over the ACDM prediction. 



4 THE 2MASS ANGULAR POWER 
SPECTRUM 

T he large local hole in the APM survey area, as suggested 
bv lBusswell et alJ <l2004l) (a «25 per cent deficiency to z=0.1 
over «4000 deg 2 ) and the K s <12.5 2MASS number counts, 
appear to imply the presence of excess power at large scales 
over the ACDM prediction. In order to investigate this possi- 
ble presence of excess power at large scales, we compute the 
angular power spectrum for \b\ ^20° 2MASS galaxies and 
2MASS N-body mocks constructed from the ACDM Hub- 
ble Volume mock catalogue. An examination of the implied 
cosmological parameters and of various possible sources of 
systematic error as well as a more de t ailed d escription of the 
method can be found in iFrith et al.l (l2004b|) . 

4.1 Method 

Following the usual meth o d (e.g. iPeeblesI Il973l : 
iPeebles fc Hauserl Il973t Il974l : IScharf et al.l Il992h . the 

angular power is estimated through a spherical harmonic 
expansion of the surface density of galaxies. The coefficients 
of this expansion are determined over the observed solid 
angle Q obs : 

where N=N ga i/Q. b a is the observed number of galaxies per 
steradian. The angular power is then determined as the ratio 
of the observed signal to the expected Poisson fluctuation, 
A/: 

Cl = A/(W^)(2/ + l)5> r|2 (5) 

m 

such that Cj=l corresponds to a random distrib ution for all- 
sky cov erage. This is described in more detail in lFrith et alJ 
<2004tft . 

4.2 Results 

The 2MASS angular power spectrum for r;120,000 K s < 
12.5, \b\ ^20° galaxies is presented in Fig. We also show 
the mean and lcr spread determined from 27 mock 2MASS 
catalogues described in section 2.4. The linear model cor- 
responding to the ACDM Hubble Volume mock catalogue 
input parameters, corrected for the window, is indicated by 
the dotted line. 

The 2MASS angular power spectrum is in reasonable 
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Figure 4. The \b\ ^20° 2MASS angular power spectrum for 
124,264 K s <12.5 galaxies. The crosses indicate the 2MASS dat- 
apoints with the shaded region and solid line indicating the la 
spread and mean power spectrum of 27 mock 2MASS catalogues 
constructed from the ACDM Hubble Volume mock catalogue. The 
model corresponding to the Hubble Volume mock catalogue input 
parameters of f2 m =0.3, f2;,=0.04, /i=0.7 and trg=0.9 is indicated 
by the dotted line. The upper panel shows the power spectrum 
normalised as in equation [5] with the dashed line indicating the 
expected Poisson fluctuation for the solid angle used. In the lower 
panel we show the 2MASS, mock and model power spectra ex- 
pressed as the fractional deviation from the mean mock power 
spectrum. 



agreement with the mock 2MASS angular power spectra al- 
though the 2MASS angular power spectrum slope is steeper 
and there is some discrepancy with the mock catalogue in 
an unbiased scenario. Therefore it appears that either there 
is an excess of power in the 2MASS catalogue over the 
ACDM Hubble Volume or there exists a scale-dependent 
bias within the scales shown. Clearly, the issue of bias is 
critical in determining the level of disagreement at large 
scales. T aking a reasonable v alue of the 7f s -band bias of 
&k = 1T (Mail er et alj 1200 3a - ) . the disagreement at large 
scales (I ^30; this corresponds to r~> 30 h~ l Mpc at the mean 
depth of the K 3 <12.5 sample) is at the «3er level. 



4.3 The relevance for a large local hole 

Using equation^ we can determine whether the possible ex- 
cess of power observed at large scales in the 2MASS angular 
power spectrum over the ACDM prediction can account for 
a large local hole in the APM survey area, via a transform of 
the angular power spectrum to the angular correl ation func- 
tion. For this we use a Bessel function transform jEfstathioul 
1990): 



Figure 5. The K s <12.5, |fe| ^20° 2MASS (crosses) and mock 
2MASS (solid line) 2-point angular correlation functions deter- 
mined via a Bessel function transform of the K s <12.5 angular 
power spectrum for 1^ I ^100. We have included a bias of bj{ =1.1 
for the mock 2MASS catalogues. The solid lines indicate power 
law fits of A=0.28, l-7=-0.71 for 6 <4.5°, A=4.4, l-7=-2.5 for 
6 ^ 4.5° where cu = A6 1 >. The small scale fit is taken directly 
from lMaller et alj i2003al) ; at large scales the fit is determined for 

e ^io°. 



ICiMW) (6) 

i 

In Fig. |SJ we present the transformation of the K s < 
12.5 2MASS and ACDM mock 2MASS power spectra shown 
in Fig. 0] to the angular correlation function, together with 
the best fit power laws used in section 3.3. Both the 2MASS 
and mock 2MASS a ngular correla t ion fun ctions are in good 
agreement with the iMaller et alj <l2003al) best fit slope at 
~ 1°; at smaller scales, the transformation becomes unreli- 
able due to the lack of angular power spectrum information 
beyond 1=100. 

The 2MASS angular correlation function is in good 
agreement with the mock 2MASS catalogues at large scales, 
although there is a small difference in slope. Computing the 
significance as in Equation [3] using a new best fit power law 
to the 2MASS angular correlation at large scales (9 10°) 
does slightly reduce the significance; the significance esti- 
mates in Tableware reduced by <0.5a in each case. 

Since the significance estimates detailed in Tableware 
only marginally effected by this change in the large-scale 
slope of the angular correlation function, it appears that 
the discrepancy found between the 2MASS angular power 
spectrum and the ACDM mock 2MASS power spectra at 
large scales is not able to account for the large local hole in 
the galaxy distribution described in section 3. The only way 
in which this local hole and the lack of excess large-scale 
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Figure 6. Count Probability Distribution Functions (CPDF) for N=0, 2, 6 and 20 for K s <12.5 2MASS galaxies (crosses). The mean 
CPDFs (solid line) and la spread from the 27 mock 2MASS catalogues are also shown, for unbiased (lower panels) and biased (upper 
panels) particles. The errorbars are Poissonian. 




Figure 7. Density Probability Distribution Functions (DPDF) for densities of p=50±5 per cent, 75±5 per cent, 100±5 per cent and 
150±5 per cent for K s <12.5 2MASS galaxies (crosses). The mean DPDFs (solid line) and lc spread from the 27 mock 2MASS catalogues 
are also shown, for unbiased (lower panels) and biased (upper panels) particles. The errorbars are Poissonian. 
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power over the ACDM prediction can be reconciled is if the 
whole local sample is significantly biased by the local galaxy 
distribution in the APM survey area and therefore that the 
2MASS catalogue is not a fair sample of the Universe. This 
is further discussed in section 6. 



5 COUNTS IN CELLS 

Our third technique to investigate the local large-scale struc- 
ture is a counts in cells analysis of the K s < 12.5 2MASS 
sample. At one level this is simply a check of the significance 
calculation used in section 3.3. However, it is also interesting 
to investigate how higher order moments in the galaxy clus- 
tering may effect the large-scale structure, and also whether 
the Hubble Volume simulation is able to reproduce the ob- 
served structure over cosmologically significant volumes. 

5.1 Method 

We sample the 2MASS survey area with a large number 
of randomly placed circular cells with an angular radius 
r. The Count Probability Distribution Function (CPDF) is 
defined as the probability of finding an exact number of 
gala xies, N, in a particu lar cell as a function of the cell size 
(e.g. lCroton et al.l20o3) . Here, we determine the CPDFs for 
N—0, 2, 6 and 20 sampling scales of < 3°. We therefore limit 
the area over which the cells are placed to \b\ ^20° and the 
2MASS K s <12.5 sample to |6| ^10° to remove spurious 
edge effects. This magnitude limit is chosen in order that 
the peak in the selection function lies within the redshift 
range of interest; at K s <12.5 the n(z) peaks at z ~0.05, 
and so this magnitude limit preferentially samples the mean 
depth of the possible local hole in the APM survey area. 

In order to probe much larger scales, of interest in this 
work, it is more useful to probe the density of galaxies rather 
than the absolute number. We define a Density Probability 
Distribution Function (DPDF) as the probability of finding 
a cell of given density (determined with respect to the mean 
2MASS density) as a function of cell size. Since we wish 
to probe the galaxy distribution at large angular scales, we 
limit the cells to |6| ^42° and the 2MASS K s <12.5 sample 
to \b\ >10° as before. For both the CPDFs and DPDFs, we 
mimic the 2MASS sample with the 27 biased and unbiased 
mock 2MASS catalogues described in section 2.3. 

5.2 Results 

In Fig. HO we show the CPDFs for iV=0, 2, 6 and 20 de- 
termined for the 2MASS and mock 2MASS samples. The 
agreement between the 2MASS and the mock galaxy sam- 
ples is excellent, although this is not unexpected since the 
CPDF is dominated by the 2-point correlation function; the 
mock galaxy sample is produced to match the observed to (9) 
at small scales. 

Moving to larger scales, Fig. |7| shows the DPDFs for 
the 2MASS and mock galaxy samples for a range of densi- 
ties, determined with respect to the mean 2MASS density. 
Again, the agreement is excellent to extremely large scales 
(#<30°). It is also interesting to note that the introduc- 
tion of bias has an extremely significant effect on the re- 
sulting DPDF. While the solid angles probed in the largest 



bin are slightly smaller than the APM survey area, the good 
agreement between the biased mocks and the 2MASS sam- 
ple confirms the significance calculation shown in Table Q 
which suggests that, when compared to the mean |6| ^20° 
2MASS number count, the observed deficiency in the APM 
survey area is not significant when compared to the ACDM 
prediction. The mock 2MASS catalogues are normalised to 
the mean number count observed by 2MASS and so any in- 
crease in the global galaxy density, as might be suggested 
by the faint Tf-band counts, increases the significance of a 
large local hole. 



6 DISCUSSION & CONCLUSIONS 

Recen t evidence jBusswell et all l2004t iFrith et alJ 120031. 
l2004af) has suggested that while the optical number counts 
over the APM survey area may be significantly less deficient 
than originally proposed, the resulting under-density might 
still present a challenge to the form of clustering predicted 
by ACDM at large scales. In this paper, we have presented 
three different methods of analysis to probe the possible ex- 
istence of this large local hole in the galaxy distribution in 
the SGC using the recently completed 2MASS survey. 

First, we determined the iG-band number counts over 
the APM survey area. In order to probe the underlying 
galaxy distribution, we compared the APM counts with a 
model guided by the Southern 2dFGRS n(z); this variable 
(f>* model provides a reasonable agreement with the corre- 
sponding isT s -band counts extracted for the Southern 2dF- 
GRS strip. The agreement between the Southern 2dFGRS 
variable 0* model and the A's-band counts over the APM 
survey area is remarkable, and suggests that the galaxy 
distribution over this ~4000 deg 2 area may be similar to 
that of the much smaller 2dFGRS Southern strip. Using the 
2MASS-2dFGRS matched n(z), this would imply a mean de- 
ficiency in the gala xy distribution of 24 per cent to z = 0.1 
with respect to the lMetcalfe et all i200ll) model used in this 
paper; taking the la lower limit best fit normalisation of 
this model to faint K-band data compiled from the litera- 
ture in the range 14< K <18 implies a mean deficiency of 
15 per cent to z=0.1. 

The issue is complicated by the b ^20° and b ^-20° 
2MASS counts which are in go od agreement w i th eac h other 
but are significantly below the lMetcalfe et all i200ll) homo- 
geneous prediction. There are three possible interpretations. 
The first is that the model normalisation is too high. How- 
ever, if the model were scaled down in order to agree with 
the mean 2MASS number density, this would require a 2.5<r 
deviation from the faint 7f-band counts; to account for the 
low APM counts entirely through a change in the normal- 
isation would require a 3.8a deviation. Secondly, the low 
cap counts might in dicate a zero-point offset between the 
iMetcalfe et alJ 1200 J) model and the 2MASS data. Invoking 
the «0.15 magnitude shift necessary to bring the cap counts 
and the model into line would comprom i se the good agree- 
ment between the 2MASS and lLovedavl (1200011 photometry 
and also between the 2dFGRS -Ks-band counts and the cor- 
responding variable (j>* models. Thirdly, the low 2MASS cap 
counts might indicate that the entire local galaxy distribu- 
tion is globally under-dense. While this might appear to be 
unlikely, the observed 2MASS counts over the APM survey 
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area suggest that large inhomogeneities exist in the galaxy 
distribution over extremely large volumes, and so perhaps 
only a few such features are necessary to bias the entire local 
sample. It is therefore not possible to rule out this final pos- 
sibility without further analysis. However, since the b ^20° 
and b ^-20° counts are similar, this position requires the 
coincidence th at we are posit ioned in the centre of this local 
underdensity jLoved av 2004). 

In order to determine the significance of the observed 
A's-band counts over the APM survey area, we calculated 
the expected variance over «4000 deg 2 considering a ACDM 
form of the 2-point angular correlation function at large 
scal es. The observed defic iency is calculated with respect to 
the lMetcalfe et all l)200l|) model, the best fit and la lower 
limit best fit normalisations of this model described in sec- 
tion 3.1, and the mean 2MASS \b\ ^ 20° number density 
(n 2mass)- In the first three cases, the observed counts rep- 
resent an ext remely rare fluctuatio n from that expected by 
ACDM. If the lMetcalfe et al.1 <200ll) model is effectively low- 
ered by 2.5cr with respect to the faint A"-band data to the 
mean 2MASS number density, then the observed counts be- 
gin to become more consistent with a ACDM form of the cor- 
relation function. Therefore, unless the faint if-band data 
are systematically too high or the galaxy distribution in the 
SGC is an extremely rare fluctuation in the galaxy density, 
then the AVband counts over the APM survey area appear 
to imply an excess of power at large scales over the ACDM 
prediction. 

Our second technique is therefore to examine the large- 
scale power in the 2MASS catalogue through a determina- 
tion of the galaxy angular power spectrum. We compare this 
to a ACDM prediction determined from the Hubble Volume 
simulations. The two are in reasonable agreement although 
there is some di screpancy in the slo pes; taking a _ff s -band 
bias of 6 k = 1.1 l)Maller et al.ll2003al ) results in a » 3(7 ex- 
cess over the mean ACDM angular power spectrum at large 
scales (I <30). In order to determine the effect of this appar- 
ent excess on the significance estimates used previously, we 
transform the angular power spectrum to the angular corre- 
lation function via a Bessel function transform. The corre- 
sponding best fit to the 2MASS angular correlation function 
at large scales (8 ^ 10°) decreases the significance estimates 
by ^ 0.5(7 compared to the ACDM angular correlation func- 
tion used previously. Therefore while there appears to be an 
excess of power at large scales in the 2MASS catalogue over 
the ACDM Hubble Volume simulation, it is not enough to 
account for the observed deficiency in the APM survey area. 
One caveat to this is that the mock 2MASS catalogues are 
constructed such that the mean galaxy density agrees with 
that of the \b\ > 20° 2MASS sample. If, as might be inferred 
from the low cap counts in Fig. 01 the local galaxy distribu- 
tion was globally under-dense with respect to the faint K- 
band counts, then the corresponding significance estimates 
would change due to the fact that the 2MASS correlation 
function is not drawn from a fair sample of the Universe. 

Our third technique is to use a 2-dimensional counts in 
cells analysis on the 2MASS catalogue and also the unbiased 
and biased mock 2MASS catalogues constructed from the 
ACDM Hubble Volume simulation (again normalised to the 
mean 2MASS galaxy density). This enables us to verify the 
significance estimates used previously, and also that the form 
of real features in the galaxy distribution are reproduced 



by the ACDM Hubble Volume simulation at large scales. 
The biased mock catalogues reproduce the observed galaxy 
distribution to very large scales (<30°). This supports the 
significance calculation determined previously (with respect 
to the mean 2MASS number density) which indicates that 
the local hole is not a challenge to ACDM if the K s < 12.5 
2MASS catalogue is a fair sample of the Universe. 

In conclusion, the issue of the large local hole in the lo- 
cal galaxy distribution has yet to be resolved. The 2MASS 
K s -b&nd number counts extracted for the «4000 deg 2 APM 
survey area indicate a clear deficiency in the local galaxy dis- 
tribution consistent with the form of the n(z) in the much 
smaller 2dFGRS Southern strip. However, to determine the 
level of the deficiency in the galaxy distribution requires an 
accurate n ormalisation of the Tf -band number count model. 
Using the lMetcalfe et ail (j200lfl model, which provides an 
excellent fit to faint if-band data compiled from the litera- 
ture in the range 14< K <18, implies a deficiency over the 
APM survey area which is at odds with ACDM and a local 
galaxy distribution which is globally under-dense. Only if 
the model is lowered by 3.8a below the faint _K"-band data 
can the normalisation account for the low counts over the 
APM survey area. The large increase in faint if -band data 
from the UK Infra-red Deep Sky Survey should help to re- 
solve this issue. 
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